Charge ordering in doped manganese oxides: lattice dynamics and magnetic structure 
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Based on the Hamiltonian of small polarons with the strong nearest neighbor repulsion, we have 
investigated the charge ordering phenomena observed in half-doped manganites Ri^A^MnOa. 
We have explored possible consequences of the charge ordering phase in the half-doped manganites. 
First, we have studied the renormalization of the sound velocity around Too, considering the acoustic 
phonons coupled to the electrons participating in the charge ordering. Second, we have found a new 
antiferromagnetic phase induced by the charge ordering, and discussed its role in connection with 
^ — , the specific CE-type antiferromagnetic structure observed in half-doped manganites. 
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Since the discovery of colossal magnetoresistance (CMR) compounds Ri_ x A x Mn03 (R = La, Pr, Nd; A = Ca, 
Ba, Sr, Pb) there have been extensive experimental and theoretical efforts to find the origin of the anomalous 
magnetotransport phenomena. Despite such intense studies, however, a complete understanding of unusual physical 
properties in perovskite manganites is still lacking. Undoped RMnC>3 is an antiferromagnetic (AFM) insulator, but 
with carrier doping of divalent A elements in R sites, the system becomes a ferromagnetic (FM) metal with the 
observed CMR for 0.2 < x < 0.5. Qualitative explanation is given by the double exchange (DE) mechanism j^]. At 
higher doping (x > 0.5), the ground state becomes again an AFM insulator A phase boundary between the 

. FM metallic and the AFM insulating ground states exists in a narrow range around x = 0.5 ||. In addition, another 
1 intriguing phase, the charge ordered (CO) state has been found to exist in insulating R^A^MnOa (RA=LaCa 
|pL PrSr Q, PrCaB, NdSr H). A direct evidence of the CO state is provided by the electron diffraction for 
La 1/2 Ca 1/2 Mn0 3 JTgT 

The CO state is characterized by the real-space ordering of Mn 3+ /Mn 4+ in mixed valent Ri^A^MnOs. The 
CO state is expected to become stable when the repulsive Coulomb interaction between carriers dominates over the 
kinetic energy of carriers [ fLl| . In this respect, the electron lattice of the CO state may be viewed as the generalized 
Wigner crystal. Furthermore, the carriers to form the CO state are believed to be manifested in some type of polarons, 
which arise from the strong electron-phonon interaction, possibly, the Jahn- Teller effect [jl2|. In fact, ordering of such 
1 1 polarons are occasionally observed in 3d transition metal oxides. The Verwey transition in Fe3 04 is a typical 
example of the real-space ordering of mixed valent Fe 2+ and Fe 3+ species. Also the formation of polaron lattice has 
been reported in La 2 _ x Sr x Ni0 4 (x = 1/3 and 1/2) @. 

Recently, Ramirez et al. || observed in Lai_ x Ca x Mn03 (0.63 < x < 0.67) that the CO transition is accompanied 
by a dramatic increase (>10%) in the sound velocity, implying a strong electron-phonon coupling. Another interesting 
aspect of the CO state is its relevance to the observed magnetic phases. In the half-doped LaCa j|,[l0| and PrSr Qj 
manganites, the CO state has been realized with FM-AFM transition. On the other hand, in PrCa [B| manganites, 
it has been found in the paramagnetic (PM) state followed by the AFM transition thereafter. However, the common 
feature is that the AFM structure of the specific CE-type || is observed in the CO state of manganites, suggesting a 
nontrivial effect of the CO state on the magnetic phase. The other noteworthy observation is the transport phenomena 
of the CO phase in the presence of the magnetic field. The high magnetic field induces the melting of the electron 
lattice of the CO phase to give rise to a huge negative magnetoresistance (MR) 

In this paper, we have studied the CO transition in half-doped manganites in the context of an order-disorder 
transition in interacting small polarons. The model seems quite plausible, in view of the fact that the transport 
behaviors in half-doped CMR compounds are qualitatively similar to the small polaron conduction in Fe304 [ |l6| . 
We have then investigated interesting physical consequences associated with the CO state. First, we have studied 
the CO effects on the lattice dynamics, taking into account the acoustic phonons coupled to the CO states. Second, 
considering the spin interaction in terms of the DE, we have studied the magnetic structure resulting from the CO 
state. It is found that the CO transition produces a large sound velocity renormalization around Tqo, and that the 
CO state favors the specific AFM structure apparently very similar to the CE-type AFM structure of half-doped 
manganites. 

To describe the charge ordering phenomena in the mixed valence systems with Mn 3+ :Mn 4+ = l:l, we start with the 
following Hamiltonian of repulsively interacting small polarons, 

n - = v E c ^ c h + *< cos 2> 5 °* Cj (1) 

<ij> <ij> 
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where < ij > denotes the nearest neighbor pair. Here, V and t represent, respectively, the nearest neighbor polaron 
interaction and the hopping parameter, both of which are renormalized by the strong electron-phonon interaction 
jl7]]. As the carrier concentration increases, the interaction between carriers will play a more decisive role in the 
system. The second term in TL C is the DE interaction with the single band, and the thermodynamic average (cos|) is 
determined by underlying spin structures. We assume that electrons are strongly trapped into small polarons so that 
there is no significant overlap of orbitals. Then the DE term can be treated as a small perturbation (V ^ t(cos|)). 
This assumption is consistent in that the ordering of mixed valent Mn 3+ /Mn 4+ atoms weakens the DE interaction, 
as the degeneracy in the ground state (Mn 3+ — O 2- — Mn 4+ and Mn 4+ — O 2- — Mn 3+ ) is lifted [Q. We also consider 
spinless electrons, since the strong on-site interaction U (~ lOeV) prohibits the doubly occupied states. Therefore the 
most important physics in our consideration is the nearest neighbor repulsion between small polarons. 

In the lowest order approximation, we can neglect the DE term in the above Hamiltonian, which then can be 
transformed into the spin one-half AFM Ising model, 

<ij> 

where we have used the pseudospin operator (cjc; = \ + \<Ji) fl9fl . The AFM transition of the above Ising model should 
be interpreted as the CO transition, and the sublattice magnetization, (er^) will correspond to the order parameter of 
a given charge order-disorder transition. The spontaneous magnetization will occur below the transition temperature 
Tn- In the mean field approximation, T/v is given by zV/A (z : coordination number) which corresponds to the CO 
transition temperature Tco- Thus below Tco, the real space ordering of Mn 3+ /Mn 4+ species takes place, minimizing 
the repulsion between localized small polarons. 

Another notable thing in the original Hamiltonian, Eq.(l) is that the CO state becomes unstable when V < i(cos|). 
Since both V and i(cos|) critically depend on the spin structure and the electron-phonon interaction of the system, 
such a situation could happen by controlling external thermodynamic variables, e.g. magnetic field or pressure. The 



subtle balance between two parameters, V and t(cos|), will actually determine the stable ground state |20|. The 
melting of the CO state induced by the magnetic field could be understood in this framework. This part will be left 
for future study. 

To investigate lattice dynamics near the CO transition, we consider the extended Hamiltonian H. cp , incorporating 
the electron-acoustic phonon coupling, 



<H> k 1 k 



where and are the momentum and the amplitude of the local vibrational mode, and the final term represents 
the electron-phonon interaction of local deformation type. We have adopted the pseudospin operator here too. By 
employing the canonical transformation, = Qj: + uo(il) the electron-phonon coupling term can be decoupled. 
Then introducing the generalized susceptibility function [jl9|, one can obtain the following renormalized phonon 



frequency, Coik) = uj (k)/ y 1 + g 2 (k)xaa(k), where x<?<?(k) is an AFM spin susceptibility of the wave vector k. 

Since we are dealing with the acoustic phonon with u>o(k) = sk, the renormalized sound velocity s in the limit of 
k — > is given by 



~S = S /v/l + 9 2 (0)Xaa(T). (4) 

Here x (T(T (T) corresponds to the uniform AFM spin susceptibility of the Ising model with the exchange interaction 
jV = j ^2j;(Vk — 2g 2 (k))e~ lk ' T (r : nearest neighbor vector). Note that, as the electron-phonon coupling is turned on, 
the AFM Neel temperature becomes reduced to Tn — zV/4(= Tco)- By using the uniform AFM spin susceptibility 
in the mean field approximation plf , the behavior of sound velocity s as a function of T can be obtained (see Fig. ^). 
Due to a cusp maximum in the AFM susceptibility, there appears a cusp minimum in the sound velocity at Tco- The 
sound velocity is softened smoothly above Tco, but becomes substantially hardened below Tco- Quite interestingly, 
the behavior of the above sound velocity around Tco is qualitatively very similar to the observation of Ramirez et 
al. H (denoted by full dots in Fig. |l|) for La!_ x Ca x Mn03 (0.63 < x < 0.67). It is noteworthy to compare the present 
sound velocity results with those in other charge ordered materials such as Lai.g7Sro.33Ni04 (Tco ~ 240K) f22|| . The 
sound velocity just near Tco i n this compound is very similar to the present calculation, but as T decreases far below 
Tco, it becomes again softened. That seems due to the residual charge hopping or other possible effects in nickelates. 
On the contrary, in the manganese oxides, the residual hopping is not available because the DE is suppressed by the 
CO. 
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The hardening of the sound velocity has been also reported for the low doping compound (La2/3Cai/3Mn03) of 
the FM metallic phase below the FM transition temperature, Tq But the hardenings of the sound velocity 

in two cases have different origins. In the low doping case, the renormalization of the sound velocity occurs due to 
T-dependent electron screening coming from the DE interaction, and so the hardening is enhanced by the external 
magnetic field which greatly influences the DE interaction p3|. On the contrary, in the high doping case, the present 
study reveals that the ordering of localized polarons is responsible for the renormalization. Accordingly, the sound 
velocity in the CO state will not be susceptible to the external magnetic field, because the DE interaction is not 
involved in the ordering and the accompanying renormalization. It is thus expected that the behavior of the sound 
velocity under the magnetic field would be very different between two cases. Different behaviors for the above two 
cases were indeed observed in the experiment of Ramirez et al. |J. Meanwhile, the stronger magnetic field will 
destabilize the CO state and definitely modify the behavior of the sound velocity in the CO state, which needs further 
experimental investigations. 

Now let's examine the magnetic structure in the CO state. Collecting relevant spin interactions, we have 

n m = y E n * n i + 1 E cos -f- c ta + Js E Si ' Sj ' ; (5) 

<ij> <ij> ~ <ij> 

where the last term is the AFM superexchange between local spins. Although the DE interaction does not contribute 
to realizing the CO state, it is the lowest order in the local spin interaction which yields a direct connection between 
the spin and the CO ordering. One can get the coupling between e g electrons and ti g local spins by expressing the 

DE, Hue — ^ J2<ij> fl + S 2 gf 3 ^ c\cj. This approximation will be valid when T is not too close to K Then 

by using the mean field approximation; Hde — 2 j 2S i S<y> ( c ! c j')^i ' one can study effects of electron orderings 
on the local spins pfj. 

Defining J' = 2 J 2S2 ( c i Ci +^ = ^ a x,iay,±az) as the spin exchange coming from Hde, we investigate how 
differently the exchange J' behaves between the disordered (uniform) phase and the ordered phase (CO). The exchange 
J{ in the disordered phase (T > Tco) corresponds to 2 J 2S2 jj J2k( c \ c k) elk T (^ : number of sites), which is nothing 
but the conventional FM interaction coming from the DE. The situation, however, becomes quite distinct in the 
ordered phase (T < Tco)- Below Tco, the charge modulation starts to occur as a result of the CO, and so one can 
express (cjci+f) as 

k QEA k 

where the second term results from the CO state, and (ctcj +( g) is related to the CO order parameter with the charge 

modulation vector Q. If the carriers in the CO phase are located in every other sites, the vector Q belongs to the set 
of A = { — (±1, ±1, ±1)}. As a result, the exchange J' is given by J' = J[ + J 2 with 

(7) 

for a given f (e l< ^' T = —1) and the lattice site vector Ri = a(lx + my + nz). Note that J' 2 gives the AFM exchange 
interaction. The strength of J 2 in the CO phase is expected to be comparable to those of J[ and J s , because J 2 is 
somehow related to the order parameter (<jj) which increases rapidly below Tco- The phase factor e l< ^' R * (— e J7r ('+ m +™)) 
in the spin exchange J 2 of Eq.(7) is of particular interest, because it gives rise to an unusual AFM spin structure. 
Even though the explicit value of J 2 is hard to evaluate, the schematic spin structure arising from J 2 can be obtained 
as in Fig. |2|(a) which shows the structure in a-b plane. It is seen that the spin structure is rather complicated with 
zig-zag shaped ferromagnetic chains which are mutually coupled in an antiferromagnetic way. Most remarkable is that 
this spin structure is very similar to the real spin structure of CE-type AFM (see Fig. ||(b)) observed in half-doped 
manganites This comparison illustrates that the AFM exchange (J 2 ) induced by the CO may be essential in 
realizing the CE-type spin structure through the competition with other exchanges such as the FM-DE (J{), and the 
already existing AFM superexchange (J s ^10 meV). 

The above arguments reflect that the DE interaction still plays a essential role in the magnetism of the highly 
doped manganites, and in this manner, one can understand consistently magnetic properties of both metallic FM 
states (0.2 < x < 0.5) and the AFM ordered electron lattices (x > 0.5) within the same framework. Finally it should 
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be pointed out that the CO transition is likely to be accompanied by the orbital ordering and the lattice distortion 
|p!8| . The latter phenomenon can be described in the present scheme by following the procedure by Bari 27 1. The 
former, however, is not considered in the present model, since the non-degenerate single band is assumed for the 
DE interaction. Nonetheless, the orbital ordering may play an important role in determining the observed magnetic 
structure [ fi8|j2§| ]. Generalized model including both the charge ordering and the orbital ordering interaction will 
certainly provide more realistic magnetic structure. 

To summarize, we have described the CO phase in R^A^MnOs in the context of the order-disorder transition of 
repulsively interacting small polarons. We have found that the drastic renormalization of the sound velocity in the CO 
transition arises from the strong coupling of acoustic phonons to the CO states. Further, we have found interesting 
implication of the effect of the CO state on the magnetic structure, which may be essential for the observed specific 
CE-type AFM structure in the CO phase of half-doped manganites. 
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FIG. 1. The behavior of the renormalized sound velocity s around the charge ordering (x=0.5). In the inset, the AFM spin 
susceptibility X<ra{T) of the Ising model is plotted. To help understanding, the comparison of present results with Ramirez 
et al.'s (x = 0.63) (full dots) are also provided (Ref. [5]). In this plot, we used g 2 (0) = 0.8zV for the parameter of the 
electron-phonon coupling strength. 
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FIG. 2. (a) The AFM spin structure (a-b plane) induced by the additional exchange J' 2 in the CO phase. White and grey 
circles denote Mn 3+ and Mn 4+ , respectively, (b) The real spin ground state (CE-type) observed in the half-doped manganites 
(Ref.[3]). In both cases, the spin structures are given by zig-zag shaped ferromagnetic chains coupled antiferromagnetically. 
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